A new exactly solvable quantum model in N dimensions 
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^ Abstract 

O^ The iV-dimensional position-dependent mass Hamiltonian 

CM fj^ -^' v2 ^W 

> 2(1 + Aq2) ^2(1 + Aq2) 

^T^ is shown to be exactly solvable for any real positive value of the parameter A. 

y—{ Algebraically, this Hamiltonian can be thought of as a new maximally super- 



r- 



integrable A-deformation of the A/^-dimensional isotropic oscillator and, from a 



^^ geometric viewpoint, this system is just the intrinsic oscillator potential on an 

^—^ A^-dimensional hyperbolic space with nonconstant curvature. The spectrum 

• • of this model is shown to be hydrogenlike, and their eigenvalues and eigen- 

.'^ functions are explicitly obtained by deforming appropriately the symmetry 

/\ properties of the A^-dimensional harmonic oscillator. A further generalization 

3 of this construction giving rise to new exactly solvable models is envisaged. 
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1 Introduction 

The iV-dimensional (iVD) classical Hamiltonian given by 

^(q, P) = r(q, p) + W(q) = ^^^^ + ^^^^y (1) 

with real parameters A > and u > and where q, p G M^ are conjugate coordi- 
nates and momenta, was proven in [1] to be maximally superintegrable. This means 
that "H is endowed with the maximum possible number of (2A^ — 1) functionally 
independent constants of motion. Explicitly, the integrals of the motion for "H are 
the ones that encode the radial symmetry of the system, namely, 

C(™)= 5]te-g,p.)^ C(„)= Yl te-9^P0', m = 2,...,N; (2) 

l<i<j<m N-m<i<j<N 

together with the functions 

I.=p^- {2Xn{q,p) - uj')ql t = l,...,N. (3) 

Moreover, each of the three sets {?/,C(™)}, {?/,C(„)} (m = 2, . . . , A^) and {li} 
{i = 1,...,N) is formed by N functionally independent functions in involution, 
and the set {"H, C^^\ C{rn)i h} for m = 2, . . . , N with a fixed index i provides the 
set of (2A^ — 1) functionally independent functions. Alternatively to the approach 
performed in [1], both sets of integrals (2) and (3) can also be obtained through a 
Stackel transform (or coupling constant metamorphosis) [2, 3, 4, 5, 6] from the free 
Euclidean motion which has been achieved in [7]. 

Maximally superintegrable Hamiltonians in A^ dimensions are quite scarce [8] 
even on the Euclidean space, and the two representative examples of this class of 
systems with periodic bounded trajectories are the Kepler system and the isotropic 
harmonic oscillator. In fact, "H (1) can be interpreted as a genuine (maximally super- 
integrable) A-deformation of the A'D Euclidean isotropic oscillator with frequency 
u, since the limit A — > of (1) yields 

n / -'- 2 I -*- 2 2 

^0=2? + 2^ ^ • 
Moreover, an important property of "H is that it can be written as 

1 ^ 

and, again, the limit A — ?■ of the integrals Jj transforms this equation into the 
separability property for the A^D harmonic oscillator. 

From a geometric perspective, T can be interpreted as the kinetic energy defining 
the geodesic motion of a particle with unit mass on a conformally flat space which 
is the complete Riemannian manifold M^, with metric 

ds^ = (1 + Aq2)dq2, (4) 



and nonconstant scalar curvature given by 

{N - 1){N{2 + 3\q^) - 6Xq^) 



R=-X 



;i + Aq2)3 



The scalar curvature R{r) = -R(|q|) is always a negative increasing function such 
that limj._^oo -R = and it has a minimum at the origin 

R{0) = -2XN{N-l), 

which is exactly the scalar curvature of the A^D hyperbolic space with negative con- 
stant sectional curvature equal to — 2A. In fact, such a curved space is the ND 
spherically symmetric generalization of the Darboux surface of type III [9, 10], which 
was constructed in [8, 11]. On the other hand, the central potential U was proven 
in [1, 8] to be an "intrinsic" oscillator potential on that Darboux space (we remark 
that in [1] we considered the Hamiltonian H with "H = kH/2 and A = l//t). 

For A^ = 3, the potential U appears in the context of the so called Bertrand space- 
times [12, 13]. These spaces are certain (3 + 1)D static and spherically symmetric 
Lorentzian spacetimes for which their bounded geodesic motions are all periodic 
and possess stable circular orbits. Therefore these spaces provide the natural arena 
for the generalization of the classical Bertrand's theorem [14] to relativistic spaces 
of nonconstant curvature, and the maximal superintegrability of the associated 3D 
Bertrand Hamiltonians has been recently proven in [15]. Again for N = 3, Ti arises 
as a particular case of the generalizations of the MIC-Kepler and Taub-NUT sys- 
tems constructed in [16, 17] (in particular, "H can be recovered by setting u = 1/2, 
a = 1 and b = X). 

Alternatively, "H can be interpreted as a position-dependent mass system in which 
the conformal factor of the metric (4) is identified with the variable mass function 

m(q) = 1 + Aq^. 

The construction and analysis of the wide range of applications of several position- 
dependent mass Schrodinger Hamiltonians can be found in, for instance, [18]-[31] 
(see references therein). In particular, one-dimensional models containing quadratic 
mass functions are considered in [25, 26] for certain semiconductor heterostructures. 

In this letter we solve the quantization problem for "H when A > 0. In the 
next section we introduce the radial effective potential associated to this system. In 
section 3 we perform the quantization by imposing the existence of the quantum 
analog of the full set of integrals of the motion. The explicit solution of the spectral 
problem is presented in section 4. Furthermore, a generalization of this approach in 
order to obtain new exactly solvable quantum models with position-dependent mass 
is sketched in the last section. 



2 The classical effective potential 

The Hamiltonian "H (1) can also be expressed in terms of hyperspherical coordinates 
r,6j, and canonical momenta pr,pe , {j = 1,...,N — 1). The N hyperspherical 
coordinates are given by the radial one r = |q| G M^ and A^ — 1 angles 6j G [0, 27r). 
In terms of them we have that 

i-i iv-i 

Qj = r cos 9j TT sin 9k, I < j < N, Qn = r TT sin 6k, 

k=l fc=l 

where hereafter any product YIT such that / > m is assumed to be equal to 1. The 
metric (4) now reads 

ds^ = {l + \r^){dr^ + rMn^), (5) 

where dQ"^ is the metric on the unit (A^ — 1)D sphere 

N-l i-1 

i=i k=i 
Thus the Hamiltonian (1) becomes 

n(r,p,.) = T(r,p.)^m = 1^^^ + ^^^^^. (6) 

where the total angular momentum is given by 

N-l j-1 

^ ' f^ sm 9k 

j=i k=i 

In these coordinates the integrals of motion C(m,) (2) adopt a compact form (the 
remaining C*^™-* and /j have more cumbersome expressions): 

N-l j-1 

Cim)= Y. pI n —2^' rn = 2,...,N- 

j=N-m+l k=N-m+l '' 

SO that C(Ar) = L^, which is just the second-order Casimir of the 5o(A^)-symmetry 
algebra of any central potential. 

The nonlinear radial potential U{r) is always a positive growing function for 
A > (see figure 1) and such that 

W(0) = 0, hm W(r) = ^. (7) 

However, we have to consider the variable mass contribution in the whole dynamics 
and this can be better understood by introducing the radial effective potential. If 
we apply the canonical transformation defined by 



Pr nf \ ^ /^^T^ , arcsinh(VAr) 



Vl + Ar2' ^' ' 2 2y/X 
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Figure 1: The nonlinear oscillator potential for A = {0, 0.02, 0.04, 0.06, 0.1} with 
cj = 1 is plotted. The upper dashed line corresponds to the isotropic harmonic oscil- 
lator with A = 0, and the limit r — )• oo gives {+oo, 25, 12.5, 8.33, 5}, respectively. 



on the ID radial Hamiltonian (6), we get 

where cat > is the value of the integral of motion corresponding to the total 
angular momentum C(jq) = L^. Hence the radial motion for the classical system can 
be described as a particle moving on a ID flat space under an effective potential 
WeflF(Q('^))- The analysis of Weg shows that, for any value of A, the energy of the 
system is always positive and it does have a minimum located at fmin, 

r^in = 2 ' ^es{Q{rmm)) = -ACat + ^ X^C% + UJ^Cn, (9) 

which for the harmonic oscillator (A = 0) gives r^j^ = y-^^ UesiQir^in)) = uj^/cn. 



Now, if A > and cat 7^ 0, then both r, Q G [0, cxd) and the effective potential 
has two representative limits: 



UJ^ 



limWeff(Q(r)) = +00, lim U,s{Q{r)) = —, (10) 

and the latter is just (7). Thus, this effective potential turns out to be hydrogen- 
like. The values for the minimum of the effective potential, r^^^ and Ues{Q{r^i-a)) 
(9) are, respectively, greater and smaller than those corresponding to the harmonic 
oscillator (A = 0) case. This effective potential is shown in figure 2. 
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Figure 2: Plot for the effective nonlinear oscillator potential (9) for A = 0.02, ctv = 
100 and u = 1. The minimum of the potential is located at rmin = 3.49 with 
Weff(rmin) = 8.2 and Weff(oo) = 25. The dashed line corresponds to the effective 
potential of the harmonic oscillator with A = with minimum Ucs{r^i^) = 10 at 
r^in = 3.16. 

3 A maximally superintegrable quantization 

In order to obtain the quantum analogue of the kinetic energy term (1) we have 
to deal with the unavoidable ordering problems in the canonical quantization pro- 
cess that come from the nonzero curvature of the underlying space (see, e.g. [6]) 
or, equivalently, from its alternative interpretation as a position-dependent mass 
Hamiltonian [18]-[31]. 

A detailed analysis of the different possible quantization prescriptions and a proof 
of their equivalence through gauge transformations will be presented in [32] . One of 
these prescriptions consists in the so called "Schrodinger quantization" , which entails 
that the quantum Hamiltonian i-L keeps the maximal superintegrability property and 
is therefore endowed with 2 A^ — 1 algebraically independent operators that commute 
with "H. This prescription allows us to use the full symmetry machinery coming from 
the harmonic oscillator, and gives rise to the following main result. 

Theorem 1. Let l-i he the quantum Hamiltonian given by 

^-2(l + Aq2)P + 2(1 + Aq2)' ^^^^ 

such that [qiipj] = ihSij. For any real value of X it is verified that 
(i) "H commutes with the following observables, 

(;M= J2(^q^P^_q.p.f^ C^^^= J2 iq.p,-q,Pif, m = 2,...,N; (12) 

l<i<j<m N-m<i<j<N 

6 



/, = p2 _ 2Ag27i(q, p) + um i = l,...,N- (13) 

where C^^^ = C{n) o-nd V. = ^ 12i=i ^i- 

(ii) Each of the three sets {U.C^'^^}, {H,C(m)} (m = 2,...,N) and {li} (i = 
1, . . . , N) is formed by N algebraically independent commuting observables. 

(Hi) The set {'H, C^™-\ C(m), h} for m = 2, . . . , N with a fixed index i is formed by 
2N — 1 algebraically independent observables. 

(iv) H is formally self-adjoint on the Hilbert space L^(]R^, (1 + Aq^)dq), endowed 
with the scalar product 



{^mx= / ^(q)<l>(q)(l + Aq^)dq. 



The proof of this result can be obtained through direct computation. Further 
details can be found in [32]. Next, under the usual differential representation given 

by 

d 

Qi = Qi, Pi = -i^<9i = -i^T^, V = ((9i, . . . , Sat), 

OQi 

the Hamiltonian (11) leads to the following Schrodinger equation 



,2(1 + Aq2) 2(1 + Aq2 

and the change to hyperspherical coordinates yields the equation 



2(1 + Ar2) \ ^"^'^ r '' ^ r^ 



h^dl - ^ '- dr + ^+ uj^r^ ^(r, B) = E^{r, 9), (15) 



where 6 = {6i, . . . , 6n-i) and L^ = C(n) is the total quantum angular momentum. 



4 Spectrum and eigenfunctions 

From the effective potential (8) one should expect that the quantum Hamiltonian 
(11) would have both a discrete and a continuous spectrum, and this is indeed the 
case. 

By taking into account that "H can be defined in terms of the first integrals 
Ii through 1-i = \ X]i=i -^i) "^^ fi'^*^ ^hat the discrete spectrum for the Schrodinger 
equation, ?^\E'(q) = i?\E'(q) (14) can be fully determined by following exactly the 
same steps as in the flat isotropic harmonic oscillator. Let us consider a factorized 
wave function together with the eigenvalue equations for the operators Jj: 



TV ^ 

^(q) = n^^(^^)' ;^/i^ = /^^^, ^ = i,...,iv. 



j=l 



Since k = Ii{qi,pi,'H), we get 

lI^M<l^) = \ {-fi'd^ + i^' - ^^E)q^) ^{q^) = f^^^{q^)■ (16) 

Therefore if we define the frequency 

n{E) = Vco^ - 2\E whenever cu^ > 2\E, (17) 

then (16) can be expressed, in a formal manner, as the Schrodinger equation of the 
one-particle harmonic oscillator 

1 ^-n'd^ + n\^) ^(q,) = ^,^{q,). (18) 

Since necessarily ip{qi) G L^(]R, dg^), the eigenvalue /ij and the wave function ip{qi) 
turn out to be 

lii = iJ,i{E,ni) = hD. ini + -j , n^ = 0,1,2. . . 

Mq^) = i^nAE,q^) = A„, (^) exp{- ^'q^ /2}H^X^q,) , /3 = ^^ (19) 

where HnX^Qi) is the nj-th Hermite polynomial and A„^ is a normalization constant. 
Now from 



/ JV \ N 

\i=l / i=l 



we obtain that 

N 



E = Mlln+—\ = hVuj'^ -2\E (n+ —) , n = ^ni, n = 0,1,2... (20) 



1=1 
which gives an explicit expression for the energies 



E = En = -ffX (n + y ] +h(n + ^\ Jh^X^ ('' + y) + '^'- ^^^^ 

Note that the degeneracy of this spectrum is exactly the same as in the A^D isotropic 
oscillator, a feature that is again a signature of the maximal superintegrability of 
the model. 

Concerning the continuous spectrum of ?^, it can be rigorously proven [32] that 
when A > 0, it is given by [|^,oo). Moreover, there are no embedded eigenval- 
ues and the singular spectrum is empty. The explicit form for the wave functions 
connected with the continuous spectrum can also be given in hyperspherical coor- 
dinates: it turns out that their radial factor can be expressed in terms of confluent 
hypergeometric functions [32]. 

On the other hand, it can also be proven that H has an inflnite number of 

2 

eigenvalues, all of which are contained in (0, ^) and their only accumulation point 



,2 



is ^, that is, the bottom of the continuous spectrum. 



4.1 Solution in hyperspherical coordinates 

Let us now consider the Schrodinger equation (15) expressed in hyperspherical vari- 
ables. This can be solved by factorizing the wave function in the radial and angular 
components and by considering the separability provided by the first integrals C(m) 
with eigenvalue equations given by 

^(r, e) = 0(r)F(^), q^)^ = c™^, m = 2,...,N. 

From it, we obtain that Y{6) solves completely the angular part and is written, as 
expected, as the hyperspherical harmonics 

pY{e) = hH{i + N- 2)Y{e), / = 0, 1, 2 . . . 

so that / is the quantum number of the angular momentum. The eigenvalues of 
the operators C(m) are related with the A^ — 1 quantum numbers of the angular 
observables as 

Ck ^ k-i, k = 2,...,N -1, cn^I, 

that is. 

Hence the radial Schrodinger equation reduces to 



-^M £. + ^^^f: - ' :: ' + ^'r^ M = m-l (22) 



U ^2(d' I A^-ld l{l + N -2) ^ , ^2.2 
2 \ ydr^ r dr r^ 

where we have introduced the frequency (17). Since this is formally the same radial 
equation as for the ND isotropic oscillator, we directly obtain the energy spectrum 
and the radial wave function: 

E = Ek,i = hnl2k + l + —], A; = 0,1,2... 

0(r) = 0fc,(E,r) = i?,A-^4 ^^(/3V), /3=a/^, (23) 



n' 



(1+^) 



where L^ ^ are the associated generalized Laguerre polynomials and B^^i is a 
normalization constant. Finally, if we introduce the principal quantum number 
n = 2k + I, that is, Ek^i = En (20) we recover the energy spectrum (21). 



5 Generalization 

A quite simple but very important remark is worth to be commented. Let us consider 
any exactly solvable constant-mass Schrodinger problem with Hamiltonian 

?)2 , ,2 

ff=^V= + yf(q) 



such that all the eigenvalues En and eigenfunctions ipn of the spectral problem 
Hipn = En ipn cau be obtained analytically. 

Then, it turns out that the exact solvability of the position-dependent mass 
Hamiltonian Hx defined by 

-^' ^2 w'f/(q) 



2(1 + A [/(q)) 2(l + Af/(q))' 

is deeply related to that of H, since the Schrodinger problem for Hx reads 

-^' V72.,.A , W'f/(q) ^,A c^A.a 



2(l + Af/(q)) ^" 2(l + Af/(q))^" " ^"' 

and this equation can be written as 

which is just the spectral problem for H with the new energy- dependent frequency 
Qx given by 

nx = V^^-2\E}i. 

Obviously, this spectral problem will present different features depending on the 
values of Qx, but in any case the exact solvability of H provides relevant information 
in order to get the eigenvalues and eigenfunctions for Hx- We stress that the addition 
of a nonzero constant to the potential f/(q) (here scaled to 1) is essential in this 
procedure. 

Therefore, the position-dependent mass system presented in this paper is the re- 
sult of taking this A-deformation approach when U{q) is just the harmonic oscillator 
potential, thus explaining the maximal superintegrability of the system. Obviously, 
new A^D radial models based on other well-known ID exactly solvable Hamiltonians 
can be constructed and will be presented elsewhere. 
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